Dynamical instability and loss of p-band bosons in optical lattices 



O 

(N 



(N 



C/2 ■ 

ctf ■ 

^— > ■ 
C ■ 
ctf ■ 
=3 



I 

c 

o 
o 



> 

in 
in 

ON 

o 
o 



X 



J. -P. Martikainen^B 
1 NORDITA, 106 91 Stockholm, Sweden 
(Dated: January 24, 2011) 

We study how the bosonic atoms on the excited p-band of an optical lattice are coupled to the 
lowest s-band and the 2nd excited d-band. We find that in some parameter regimes the atom-atom 
interactions can cause a dynamical instability of the p-band atoms towards decay to the s and ri- 
bands. Furthermore, even when dynamical instability is not expected s- and d-bands can become 
substantially populated. We also find that, the stability properties of the excited bands can be 
improved in superlattices. 
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I. INTRODUCTION 

The study of cold atoms in optical lattices has seen a 
dramatic experimental progress in the recent past [l|, 0] ■ 
Due to realization of optical lattices, low densities, and 
low temperatures, a fantastic degree of control has been 
obtained which has made detailed studies of strongly 
correlated quantum systems possible. For example, the 
Mott-superfluid transition has been successfully ob- 
served in optical lattices. 

The early experiments were confined to the lowest 
band and while increasing interactions can make excited 
band populations non- negligible [5| , the lowest band still 
dominates. Experimentally, however, atomic population 
residing on the excited bands can be obtained by cou- 
pling atoms from the lowest band to the excited bands. 
This has been experimentally demonstrated by accelerat- 
ing the lattice for a short period Q or by coupling atoms 
from the lowest band Mott insulator into the first ex- 
cited p-band of the lattice via Raman transitions between 
bands [7] . In this latter study, it was in particular found 
that the lifetimes of p-band atoms can be considerably 
longer than the tunneling time-scale in the lattice and 
it was also possible to explore how coherence of bosons 
on the excited band was established. Very recently p- 
band bosons and their orbital ordering in the superfluid 
phase were studied experimentally in a bipartite optical 
lattice |8(. In this experiment atoms were transferred to 
the p-orbitals by deforming the superlattice in such a way 
that atoms originally in the lowest orbital of certain sites, 
could tunnel to the p-orbital of the neighboring site. 

The exciting possibilities in the physics with higher 
band atoms has naturally attracted also theoretical at- 
tention. P-band bosons have been studied in equilib- 
rium Q and as a platform to realize interesting correlated 
quantum states [!JJ-[l3| . Atoms on the p-band have also 
been found to have interesting rotational properties [l4| . 
Furthermore, there are theoretical studies where multi- 
orbital description of bosons is found necessary since it 



can induce important corrections to the simplest lowest 
band Hubbard model [HI, [l6| • 

The purpose if this article is to point out the exis- 
tence of dynamical instability which can adversely af- 
fect the stability of the p-band atoms in the superfluid 
phase. This instability is caused by the scattering of p- 
band atoms into s- and d-band atoms. While the anhar- 
monicities in an optical lattice can make this process off- 
resonant for all quasi-momenta, interplay between band 
structure and atom-atom interactions can still induce in- 
stability. We also find that even in the absence of dy- 
namical instability, the s- and <i-bands can become sub- 
stantially populated. Furthermore, dynamical instabities 
can be removed in superlattices 

The paper is organized as follows. In Sec. [TT] we out- 
line a theory which captures the essential physics of the 
loss of p-band atoms to s and <i-bands Furthermore, in 
this section we discuss realistic magnitude of various pa- 
rameters occurring in our theory. In Sec. IIIII we derive 
and solve the Bogoliubov-de Gennes equations for the 
collective modes and find the possibility of dynamical in- 
stability. In Sec. MI Al we point out the importance of 
condensation for the occurrence of dynamical instability, 
but find that s- and d-bands can be substantially pop- 
ulated even in the absence of dynamical instability. In 
Sec. IIIIBI we analyze similar processes in superlattices 
and find that there the dynamical instabilities can be re- 
moved. We end with few concluding remarks in Sec. IIVI 



II. FORMALISM 

Let us consider bosonic atoms at zero temperature in 
an optical lattice 



V(r) = V L sin 2 
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where d = A/2 is the lattice constant. To simplify the 
formalism we assume a similar lattice potential as in the 
experiment by Miiller et al. Q so that the lattice is deep 
in y- and z-directions with depth Vd — 55 Er, where 
Er = h 2 k 2 /2m is the recoil energy, m the atomic mass, 
and the wavevector k is related to the laser wavelength 
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FIG. 1: A band structure in units of Er for the three lowest 
bands for a lattice depth of Vl = 17 Er. 



(A = 843 nm) through k — 2tt/X. Throughout our cal- 
culation we will further assume that we are dealing with 
87 Rb atoms. Since the lattice is deep is two directions 
the dynamics in these directions is frozen and the sys- 
tem is effectively one-dimensional. For the atomic wave- 
functions along y- and z-directions we use the Gaussian 
ground states corresponding to the harmonic oscillators 
along y and z. Using this effectively one-dimensional 
potential together with the usual kinetic energy oper- 
ator K = -h 2 V 2 /2 m, we get an ideal Hamiltonian 
Hiat = K + V(r) from which we numerically calculate the 
dispersions [173 E s (q),E p (q), and Ed(q) for the 3-lowest 
(s, p, and d) bands as a function of quasi-momentum q. 
Typical band structure is shown in Fig. Q] 

Atoms interact through the two-body interaction [l8[ 
g/2 J dr^T(r)^T(r)^( r )^j(r), where V>(r) (V'U 1 *)) is the 
bosonic annihilation (creation) operator and g — 
4irh 2 a/m. For 87 Rb atoms the scattering length a rj 
llOao- Our main interest is to explore the role of the 
processes that transfer atoms from the condensed p-band 
to the lowest s- and the second excited d-bands. For this 
purpose we expand the field-operators in terms of the 
field-operators ip a for the three different bands. 

-0(r) = (j)p(x)ip p + (/> s (x)4> s + 0d(x)$d. (1) 

For the pure Bose-Einstein condensate in the p-band 
at quasi-momentum q, (f> p (x) corresponds to the Bloch- 
wavefunction u^(x) at quasi-momentum q. The energy 
of the p-band atom is minimized with the choice q = n/d. 
Likewise the functions </> s (x) and 0<;( x ) are related to the 
Bloch-wavefunctions u*' d (x) which dominate the dynam- 
ics at short times. Substituting the above expansion, we 



find a generic interaction term relevant to our discussion 
Hi = ^Y^liplippipp + 2n p (g sp n s + g pd n d ) (2) 
+ 2g ppsd n p 

+ 9ppsd 

$d$p$p + h - c ) ' 

where n a = ip^ipa- We further assumed that s- and ri- 
bands are almost unpopulated by dropping terms which 
were of higher order in ip s and ipd and dropped the 
far detuned processes 2p <h> 2s and 2p o 2d. The 
form of the interaction Hamiltonian here is quite gen- 
eral and only assumes the validity of the expansion 
into 3 modes and the weak population of 2 of these 
modes. The various interaction parameters are defined 
by the integrals # Q/3 = / dr|^ a (x)| 2 |^(x)| 2 and g ppsd = 
J dr(f> p (x) 2 (j) s (x)* (f)d(x)* ■ It should be noted that in 
estimating the interaction strengths, rather than using 
Bloch-wavefunctions for Q (x), we could also have used 
the better localized Wannier states, or their harmonic ap- 
proximations without affecting our results qualitatively. 
Even the quantitative differences are not dramatic since 
the estimated interaction strengths are not very different 
in the parameter regimes we consider. 

The non-interacting part of the Hamiltonian is simply 
given by Hq — E s n s + E p n p + E^rid, so that our total 
Hamiltonian is the sum H = Hq + Hj. Note that here 
we also make an approximation that the atoms are suf- 
ficiently localized so that tunneling dynamics plays only 
a minor role at the timescales of interest. However, tun- 
neling effects are to some extent still included since they 
influence the band structure and are consequently incor- 
porated into the energy levels E n . This approximation 
is worst for the d-band (which is initially unoccupied) 
whose bandwidth is by far the largest. However, tunnel- 
ing processes in the <i-band would mainly give rise to (in 
practice) irreversible atoms loss from the p-band since 
the more mobile <i-band atoms are free to relax towards 
their zero quasi-momentum minimum, where the process 
2p s + d is very far detuned. 

The detuning 5(q, k) = 2E p (q) -E s (q-k)- E d (q + k) 
which is related to the process 2p O s + d described 
by the terms •i^l^'tpp'tpp + h.c. in the Hamiltonian, plays 
an important role for the stability of the system. If one 
applies Fermi's golden rule to study the life-time of p- 
band atoms @, the lowest order processes are no longer 
energetically allowed if the lattice is deeper than about 
Vl = 18 Er, since then S(q,k) is non-vanishing for all q 
and k. If the process is always detuned for the p-band 
atoms at q = ir/d the detuning is minimized with the 
choice k — 0. Even if the detuning is zero somewhere, 
for the lattice depths considered here, the resonance is 
still located close to k = so we choose k = without 
introducing large additional uncertainties when comput- 
ing the coupling coefficients of the interaction Hamilto- 
nian. For very deep lattices the detuning approaches a 
constant value and in this sense the system differs from 
a harmonic oscillator even in this limit (T3 |. However, 
relative to the onsite trap energy scale huj the detuning 
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FIG. 2: Typical behavior for the relevant parameters as a 
function of lattice depth Vl- In (a) we show the bandwidths 
(solid line for the d-band, dashed for the p-band, and dot- 
dashed for the s-band). The bandwidth is defined as a differ- 
ence between the maximum energy within the band and the 
minimum energy within the band. In (b) we show g pp (solid), 
g ap (dashed), and g p d (dot-dashed), and (c) shows the behav- 
ior of the coupling g pps d responsible for population transfer 
away from the p-band. Finally (d) shows the minimum and 
maximum values of the detuning S(q, k). When the minimum 
value of 5(q, k) > the process 2p f-> s + d is detuned for 
all quasi-momenta. All energies in the y-axis are expressed in 
terms of the recoil energies Er. 



does approach zero since w oc \/Vl ■ 

In Fig. [5] we summarize the behavior of various param- 
eters of our model as the lattice depth is varied. It is to 
be noticed that coupling coefficients appearing in Hi are 
typically of the same order of magnitude as the detuning 
2E p (q) — E s (q) — Ed(q). For this reason the system can 
be poorly described both in the narrow band as well as 
in the wide band limit. 



III. NONLINEAR CONDENSATE DYNAMICS 

In the previous section we derived a minimal theory 
to describe how p-band atoms are coupled to the s- 
and d-bands. Let us now investigate what this model 
implies for the dynamics of pure Bose-Einstein conden- 
sates. In this limit we can describe the system in terms 
of complex amplitudes ip s (t), ipp(t) , and ipd(t) which 
are normalized to the number of particles per site i.e. 
\ip s \ 2 + \ip P \ 2 + \*Pd\ 2 = nr- This description also follows 
when the onsite wavefunctions are described by coher- 
ent states and provides an accurate picture of the system 
deep in the superfluid regime [l3[ . The equations of mo- 
tions for the order parameters are given by 



.dip s 
'~dt 



[E s + 2g sp n p ) ip a + 2g P p sd n p ipd+gpp S d'>Ppipd-> ( 3 ) 



dip 

i-Qjr = (E P + g pp n p + 2g sp n s + 2g pd n d )ip p 



and 
.dtpd 



+ 2g Pps d{i>d?Ps + c.c)ip p + 2g ppsd ip p ipd4>s, (4) 



(E d + 2g pd n p ) i-!d + ^9 P p S dn p iis+gppsdi l p ^* s - (5) 



In order to study the stability of the system we now derive 
the relevant Bogoliubov-de Gennes equations. Since the 
system is initially prepared in the p-band we have (at 
short times) ip s (t) — 5%l> s e~' lE '^ , ipd(t) = Sip^e - ^^ , and 
Tpp(t) = [V^r + Sipp] e~ lE >> t , where symbol S is used to 
indicate a small quantity. Furthermore, we defined E' s = 
E s + 2g sp n Tl E' d = E d + 2g pd n T , and E' p = E p + 2g pp n T . 
By keeping only terms which are lowest order in 6ip a we 
get a pair of equations 



iSip s = 2g ppsd n T e l{E ^ ^Sipd + g P p sd n T e l<5 '*(5Vd 



and 

i6ip d = 2g ppsd n T e +i{E '<t- E '° )t 5ij} s + gppsdn T e~ lS,t 8^* s , 

where 8' = 2E' p -E' s -E' d . The first terms on the right 
of these equations can be dropped since they are oscil- 
lating rapidly and their time average vanishes over the 
timescales of interest. By defining 



and 



5il> s {s) = [u s e~ 1 ^ + v*e +luJt ] e~ lS ' /2t (6) 



5W) = Wde- iut + v* d e +iut ] e- lS '' 2t (7) 



we find the Bogoliubov-de Gennes eigenvalue problem for 
collective modes 



cofjX = MX 



(8) 



described by their frequency u> as well as amplitudes u 
and v a . Here A T = (u Sl v s ,u d ,Vd), 
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(9) 



and 



—S'/2 gppsdriT 

M= I ° ~ 5 '' 2 9PP X% n I 111,1 

gppsdriT -d/2 1 

gppsdn-T -S'/2 



This eigenvalue problem has solutions 



1 



(11) 



Importantly, these solutions are imaginary and indicate 
dynamical instability if \8'\ < 2g pps driT- In a parameter 
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regime where dynamical instability exists a system ini- 
tially prepared on the p-band will lose atoms to the s- 
and d-bands at a rate given by 

T - ^\6' 2 -4(g ppsd n T n (12) 

It is also easy to see when this instability is more likely to 
occur. For relatively deep lattices S = 2E p — E s — E4 > 
since the anharmonicities of the lattice potential shifts 
the d-states lower most relative to the harmonic oscillator 
energy levels. On the other hand the effective detuning 
is given by 5' = S + 2nT(g pp — g sp — g P d) and will move 
closer to zero if g pp — g S p — 9pd < 0, a condition which usu- 
ally holds since all the coupling coefficients have similar 
magnitudes. Then whenever 5' becomes small relative to 
QppsdnT (which increases with lattice depth) dynamical 
instabilities can be expected. 

In a Fig. [3] we show the loss rate of the p-band atoms as 
a function of lattice depth for two different onsite atom 
numbers. This figure suggests that for small onsite atom 
numbers a dynamical instability can be present with re- 
alistic trap parameters above Vl ~ 19 Er and below 
Vl = 25 Er at which the system is expected to be in 
the superfluid regime Q. For higher atom numbers the 
region of instability is greatly increased and extends into 
the p-band Mott-insulating region where our assumption 
of coherent states is clearly invalid. Furthermore, we 
find that the timescale for the instability is substantially 
less than the timescale for the tunneling of the p-band 
atoms justifying our earlier approximation to ignore tun- 
neling dynamics. Interestingly, the dynamical instability 
is expected for lattice depths at which the lowest order 
Fermi's golden rule becomes inapplicable [9]. 

We have also solved the Eqs. ©-f5|) numerically even 
by modifying them to include interactions between s- 
and d-band atoms. We compare the Bogoliubov- de 
Gennes approach with the full Gross-Pitaevskii equations 
in Fig. 01 As is clear the agreement between the nu- 
merical and analytical Bogoliubov- dc Gennes solution 
is very good when the p-band population is dominant 
and the small difference is mainly due to the facts that 
short time dynamics has some sensitivity to phases of 
the amplitudes and that numerics also included many far 
off-resonant processes giving rise to fast oscillations on a 
much shorter timescale. These processes were ignored in 
the Bogoliubov- de Gennes approach derived earlier. If 
in the dynamically unstable region the initial state has 
a substantial s-band population the time-evolution be- 
comes more sinusoidal, but even then the maximum ri- 
band population is substantial and cannot be ignored. 

The numerical solutions also reveal a revival in the p- 
band population after the initial exponential loss. How- 
ever, we believe that this revival is unlikely to persist in 
a realistic system, due to the high mobility of the d-band 
atoms. It is more likely that those atoms coupled to the 
d-band will be irreversibly lost from the p-band. Also, in 
a trap the density distribution of the more mobile d-band 
atoms will be much broader since the length scale for the 
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FIG. 3: Loss rate Y of the p-band atoms as a function of 
the lattice depth Vl for n-r = 1 (dashed) and tit = 2 atom 
per site. For comparison, with a dot-dashed line, we show 
the result where onsite wavefunctions were assumed to be 
harmonic oscillator wavefunctions and tit = 2. The loss rate 
is given in units of 1/r = Er/H. For 87 Rb atoms r is about 
49 fis. 
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FIG. 4: Comparison between the Bogoliubov- de Gennes re- 
sult (dashed line) for the s- and d-band populations as a func- 
tion of time and the solution to the full Gross-Pitaevskii equa- 
tions (solid line) . The lattice depth was chosen as Vl = 20 Er 
and rir — 1 so that the system was predicted to be in the un- 
stable regime. Initial state had a 10 -6 population on s- and 
d-bands (to seed the dynamics) and the phases of the ampli- 
tudes were taken to be zero. 



density distribution in a harmonic trap is proportional 
to the i?y 4 , where Bd is the d-band bandwidth. For this 
reason, many d-band atoms end up spatially separated 
from the s- and p-band atoms residing closer to the trap 
minimum. 
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Role of number fluctuations 



0.3n 



In the previous section we described the system by 
assuming pure condensates in all bands and found a pos- 
sibility of dynamical instability. The question then arises 
that what role did the assumption of condensation actu- 
ally play in our result? Alternatively, deep in the Mott- 
insulating regime with tit = 2 atoms per site we can at 
short times use an ansatz 



e 0.25- 



|V(t))=V-«j(*)|l,0,l)+Vp(*)|0,2,0), (13) 

where \n s ,n p ,n d ) is the Fock state. To the lowest or- 
der in tunneling the many-body wavefunction factorizes 
into single site solutions, so the solution in a single site 
amounts to a solutions throughout the insulator. Due 
to the presence of the |1, 0, l)-state we must furthermore 
include the density-density interaction 2g sd n s n d between 
s- and d-bands into our model. In the earlier mean-field 
theory this term was a small quantity which could be 
safely ignored. We then get (at short times) 

H\i>(t)) = [(E„ +E d + 2g sd )ip sd + g ppsd V2^ |1, 0, 1) 
+ (2E P + g pp )ij} p + g ppsd V2i(j sd |0,2,0). (14) 

In this case we can easily solve the equations of motion 
for ip p {t) and ip s d(t) and find a simple Rabi-problem with 
purely real eigenenergies 

E ± (n T = 2) = 2E P + E s+E d +9 PP + 2g sd 



± 2 V {2E p ~ E s-E d + g PP - 2g sd y + 8g 2 ppsd 

and an initial state \ip(t = 0)) = |0,2,0) evolves in such 
a way that the maximum population of the |1, 0, 1) state 
is given by 



-ysd 



y pP sd 



(2E P - E s - E d + g pp - 2g sd ) 



a 2 ' 



(16) 



In Fig.[5]we show the maximum population of the |1, 0, 1) 
state which gives in indication of how reliable a pure 
p-band theory can be. Since our ansatz with a fixed 
atom number is expected to be reasonable only for deep 
lattices, in this figure we estimated the interaction pa- 
rameters by approximating the lattice site with a har- 
monic oscillator. As is clear from this figure, substan- 
tial fraction of the 87 Rb population can reside outside 
the p-band. In fact, P^/ax can even reach unity when 
2E p — E s — E d + g pp — 2g sd — which formally happens 
somewhat below the lattice depth of Vl = 20 Er. How- 
ever, in this regime stronger tunneling invalidates the 
simple three-state description as well as the assumption 
of fixed onsite atom number. 

With exactly ut — 1 atoms per site interactions do not 
contribute and instability is not expected. These sim- 
ple exercises indicate an important role of onsite number 
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FIG. 5: The maximum population of the |1,0, l)-state when 
the number of atoms tit — 2 as a function of lattice depth. 
The interaction parameters were estimated by approximating 
the lattice site by a harmonic oscillator while the energy lev- 
els were computed with a band-structure calculation with a 
quasi-momentum q — n/d. 



fluctuations in feeding the dynamical instability. What 
happens if the system is neither a Mott-insulator nor a 
pure condensate, is unclear. Presumably, the system is 
more stable than our calculations assuming a pure con- 
densate suggest, but the stability cannot be taken for 
granted. 



B. Dynamical instabilities in superlattices 

Since p-band bosons have also been studied in a bi- 
partite optical lattice [1, [l9j . let us briefly discuss the 
stability properties of p-band atoms in superlattices. We 
will again assume the same deep optical lattice along y- 
and z directions, but for the lattice along a;-direction we 
use a potential [2(| 



V(x) = Kssin 2 



V L sin 2 



which is characterized by the depths Vs and Vl of the 
short and long lattices respectively. 

In Fig. [6] we show examples of the energy levels to- 
gether with the three lowest eigenstates of this potential 
while in Fig. [7] we show an estimate of when dynamical 
instabilities might be expected. As is clear from these 
figures dynamical instabilities might occur only when 
Vl/Vs is fairly large and the eigenstates are well local- 
ized to the deep sites. When Vl/Vs becomes smaller 
instability is quickly suppressed. The reason for this is 
twofold. First, as Vl/Vs becomes smaller the energy 
levels become very different from the energy levels of a 
harmonic oscillator at deep sites and this implies an in- 
crease in the detuning 2E p — E s — E d which reduces the 



6 




-0.5 0.5 1 -1 1 

x/d x/d 




-0.5 0.5 1 -1 



x/d x/d 

FIG. 6: Energy levels in units of Er (together with the lattice 
potential) and wavefunctions in a superlattice. We fixed the 
short lattice depth to Vs = 15 Er while Vl = 1.5Vs in (a) 
and (b) or Vl = 0.8Vs in (c) and (d). In (b) and (d) the dot- 
dashed line is the wavefunction for the s-orbital, solid line for 
the p-orbital, and the dashed line for the d-orbital. 

likelihood of dynamical instabilities. Second, while for 
large Vl/Vs all the wavefunctions are well localized to 
deep sites and their overlaps are large, for smaller Vl/Vs 
the d-orbital becomes peaked in the shallow sites. When 
this occurs, the overlap between the d-orbital and the 
s- and p-orbitals is drastically reduced and g pps d is re- 
duced by more than an order of magnitude. From this we 
can conclude that supcrlatticcs can enhance the stabil- 
ity properties of p-band bosons. However, this happens 
at the cost of greater hybridization of the p-orbitals in 
deep sites with the s-orbitals in the shallow sites. In fact 
this hybridization played an important role also in the 
experiment by Wirth et al. Q . 

IV. CONCLUSIONS 

We have pointed out a dynamical instability which can 
affect the stability of p-band bosons in the broken sym- 
metry phase. All our numerical examples assumed 87 Rb 
atoms and consequently some details are expected to be 
different for different atoms. In particular, the scatter- 
ing length will be different for different atoms or it can 
be tuned with Feshbach resonance and this can have a 
dramatic effect on how big a role the coupling to the s- 
and (i-bands plays. For convenience we restricted our dis- 



cussions to an effectively one-dimensional system. If the 
dimensionality is increased we have several degenerate 
flavors in the p- and d-bands and many more interaction 
channels. However, the Gross-Pitaevskii equations anal- 
ogous to Eqs. ©-(O would still look very similar and the 
stability analysis is likely to reveal similar instabilities as 




FIG. 7: Loss rate of the p-band atoms in a superlattice due 
to the dynamical instability as a function of Vl/Vs when the 
short lattice depth was fixed to Vs = 15 Er. The coupling 
coefficients and band energies were estimated from the q = 
eigenstates and the loss rate is given in units of 1/r = Er/Ti. 



we found here for one-dimensional systems. 

We found that while p-band leaks into s- and (i-bands 
even in the Mott-insulating regime, for the dynamical 
instabilities to be present it was important that the sys- 
tem is not a Mott-insulator. However, even in a Mott- 
insulating regime for p-band atoms, these atoms are cou- 
pled to s and d-bands. Since the bandwidth for the ri- 
band atoms is dramatically larger than for lower bands, 
d-band atoms would typically be deep in the superfhiid 
regime. Under such circumstances the p-band atoms are 
coupled to a coherent d-band and this coupling will in- 
duce number fluctuations also on the p-band. How this 
effect affects the superfluid Mott-insulator transition for 
p-band atoms has not yet been explored. 

Finally it should be noted that the stability proper- 
ties of fermionic atoms are likely to be better than those 
of bosonic atoms. This is because with fermions one 
can populate the p-band by hrst filling the lowest band. 
When the lowest band is filled, the coupling between p, 
s, and d-bands becomes Pauli blocked. 
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